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The precise measurement of Doppler shift and dispersion (frequency broadening) 
of a radio echo returned from a planetary surface can yield useful information on planetary 
orbits and surface characteristics. Methods employed to measure Doppler displacements 
include banks of narrow receiving filters, phase-locked oscillator loops, and digital sampling 
and calculational techniques. The latter, while complex and sometimes slower, offer greater 
flexibility and in some cases enable measurements not otherwise possible. 



1. Introduction 

As the sensitivity of radar systems lias grown in 
recent years through higher power transmitters, 
larger antennas, and low-noise receivers, so has the 
sophistication of the assigned tasks and measurement 
techniques. In the beginning it was sufficient merely 
to establish the presence of an echo with reasonable 
certainty. As stronger signals were obtained, 
precise measurements of time delay enabled accurate 
determinations of distance to the echoing body. 
With even stronger signals it was possible to see the 
degree to which the returned power was dispersed 
in time, and thus to study the radio scattering law of 
the distant surface. 

The use of transmitters controlled by extremely 
stable frequency sources has also permitted the ex- 
ploration of the frequency shift and dispersion of the 
returned signal. This coordinate, which is in a sense 
conjugate to the time dimension associated with 
radar range, contains useful information on the state 
of motion of the echoing body. In particular, if the 
frequency dispersion occasioned by surface roughness 
and rotation is small enough, it is possible to make 
measurements of the relative orbital motions of the 
earth and distant object with an accuracy which ex- 
ceeds 0.1 m/sec (see remarks relating to figures 3 
and 4) . Such precision can significantly improve the 
accuracy with which the orbits of the planets are 
known. 

The frequency dispersion of the returned echo de- 
pends on the degree of surface roughness and the 
amount of rotation of the target. In some cases 
where the rotation is known from optical observa- 
tions, or where the signal is sufficiently strong to per- 
mit unambiguous determination of the spectral limits 
corresponding to the limbs, the surface roughness 
may be deduced solely from the shape of the spec- 
trum. Where this is not possible, an independent 
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measure of surface roughness from the observed time 
dispersion may conversely be combined with the spec- 
tral data to yield a value for rotation. 

An even more powerful method for studying the 
detailed scattering characteristics of a distant surface 
involves obtaining a frequency spectrum at each of 
a number of values of time delay. In this way it is 
possible to develop a map of surface reflectivity with 
detail which is partly determined by the frequency 
resolution available. 

It is clear, therefore, that many of the objectives 
of radar astronomy are dependent on measuring ac- 
curately the amount by which various portions of 
the returned echo power have been displaced in fre- 
quency. An attempt will be made in the following 
section to describe some of the techniques which have 
been used. 

2. Techniques 

The techniques which have been used for spectral 
analysis may in general be broken down into those 
which are primarily analog and those which depend 
heavily on digital processing. The choice of a proper 
mix between the two for a given experiment will de- 
pend upon the echoing characteristics of the target 
and upon the speed and accuracy of the components 
available. Perhaps the most obvious technique is 
simply to use one or more very narrow-band filters 
to sequentially analyze portions of the received spec- 
trum, as was done by Fricker et al. [1958], for lunar 
echoes. In this case a single 1 c/s filter was used, 
and the narrow-band tape-recorded data were rerun 
many times at different speeds until a complete spec- 
trum of the desired region was obtained. The resolu- 
tion of this system was limited by its stability to 
about 1 c/s. 

In another approach to the problem used by 
Browne et al. [1956], the autocorrelation function for 
the temporal fading of the received power was calcu- 
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lated. By taking the Fourier transform of the meas- 
ured autocorrelation function, these workers were 
able to obtain a power spectrum of the received sig- 
nal, without the stability restrictions imposed by 
direct analog filtering. Since they worked with a 
detected signal, however, their spectrum was folded 
and suffered some loss of information. 

More recently, Goldstein [1961] and Weinreb 
[1963] have cleverly adapted the autocorrelation 
approach to spectral analysis to digital manipula- 
tion. By severely limiting the amplitude excursion 
of the received signal to a point where it could be 
described as simply positive or negative (one-bit 
quantization) the calculation of the autocorrelation 
function becomes very rapid. Both workers analyzed 
bandwidths of about 35 kc/s sampled at about 70 
kc/s using simple logical circuits and delay lines to 
calculate between 20 and 45 channels of delay. 
Calculation was rapid enough to permit the data 
to be processed and integrated as rapidly as it was 
received. Necessary weighting and Fourier trans- 
formation (see eq (4)) of the calculated autocorrela- 
tion function was performed in a few minutes at 
the end of the run. 

The advantages of this novel technique of spectral 
analysis are great while the drawbacks are small. 
The advantages consist in the high degree of stability 
afforded by a digital system which in effect synthe- 
sizes a large number of identical filters. Integration 
can proceed over a number of days in order to 



improve the statistics of a noisy signal, while the 
precision with which the various analyzing filters 
can be synthesized enables one to detect very small 
differences between adjacent portions of the spec- 
trum. The drawback lies chiefly in the decreased 
final-output signal-to-noise ratio produced by the 
amplitude clipping. If the processed signal has 
Gaussian statistics, however, this loss may be made 
up completely by merely extending the observing 
period by a factor of approximately two. A block 
diagram of a typical system of this sort is shown in 
figure 1. 

The processing may be described mathematically 
as follows [Weinreb, 1963]. Let X(t) be the voltage 
waveform at some convenient (band-limited) inter- 
mediate frequency. Define: 



Y(t)=l,whenX(t)>0; 
F(0=~1, X(t)<0. 



(1) 



Sample to obtain Y(t) at intervals: At=l/2B, 
where B= bandwidth of signal. Calculate N values 
(channels) of delay for the one-bit autocorrelation 
at multiples of At: 



p' n (7iAt) ~J^y(kAt)y(kAt+nAt), 



(2) 
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Figure 1. Block diagram of the one-bit autocorrelation system. 
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where the integration in any channel has been al- 
lowed to proceed for a time T=KAt. Correct the 
resulting values of the one-bit autocorrelation 
function according to the relation of VanVleck 
[1943]: 



P n (nAt) =sin [ (tt/2) p^nAt)}; 



(3) 



the spectral estimate is then found from the Fourier 
transform of the autocorrelation: 



P(jO=2A*[po+X; PnW n cos (2ir>A*)], (4) 

n-l 



where W n is a weighting function chosen so that 
W =l and W n =0 when n>N-l. The value of N 
and the values of Wn for 0<n<iV — 1 are chosen 
according to the desired frequency resolution and 
sidelobe level; generally, the greater the resolution, 
the greater will be the fluctuation in the spectral 
estimate for a given amount of observed data 
[Blackman and Tukey, 1958]. The resulting spectral 
estimate will be normalized; it is a consequence of 
the one-bit quantization that absolute signal level 
is not preserved. 

A somewhat different approach has been used by 
Pettengill and Henry [I960, 1962] to study the 
moon. In this method a series of unnormalized 
spectra is desired of a nongaussian signal. The 
spectra have been calculated directly from digitalized 
samples of two waveforms representing the sine and 
cosine components of the returned echo signal. 



Amplitude clipping was negligible since 9-bit quanti- 
zation (512 levels) was used to sample each com- 
ponent. In this measurement the incoming signal 
represented the echo of a short radar pulse, rather 
than a continuously received waveform as in the 
earlier cases. Tims the sampling interval is set by 
the interpulse interval of the radar, which in the 
experiments described above was about 40 msec. 
Such a long sampling interval was possible because 
the returned spectrum was never wider than about 
10 c/s, while the sampling rate permitted an un- 
ambiguous measurement over a window of 25 c/s. 

Figure 2 shows a block diagram of the setup. 
Fourier transformation was carried out over 20 
individual four second segments of the received 
signal in a digital computer associated with the radar. 
The squared magnitudes of the transforms were then 
summed to yield the spectral estimate. (This corre- 
sponds to a linear taper on the W n in (4)). Because 
of the length of the digital word, it was not possible 
to carry out the spectral processing as the data were 
taken. Records were made on either paper or 
magnetic tape and reduced at a later time. The 
frequency stability of the equipment was better than 
0.01 c/s over several minutes of time, permitting 
useful resolutions as low as 0.1 c/s and a doppler 
measurement accuracy approaching 0.01 c/s. Fig- 
ures 3 and 4 show sample echo spectra taken from 
the Moon and Venus, respectively. The resolution 
is about 0.25 c/s, and the theoretical standard devia- 
tion of the fluctuation is 22 percent of the mean 
value of the estimate at each frequency. (This 
corresponds to 40 Blackman-Tukey degrees of free- 
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Figure 2. Block diagram of the 9-bit digital system used to obtain power spectra. 

1027 



(a) 

t = 2.100ms 




(b) 

t = 2.600ms 



(C) 
t=3.IOOms 





+ 865.0 +870.0 



+ 875.0 



+ 865.0 



+ 870.0 



+ 875.0 



FREQUENCY (c/s) 



Figure 3. Lunar echo power spectra at 4 intervals of range taken by Millstone Radar of MIT at 19h 00m 19s UT, June 20, 1961. 



dom.) Note that although the spectra repeat at 
multiples of the sampling frequency, the correct 
region can be selected easily from a priori knowledge. 
In the cases shown it is quite easy to locate the 
spectral center to within 0.1 c/s; since the operating- 
frequency was 440 Mc/s, the radial component of the 
relative orbital motion of the Moon or Venus can 
be determined to better than 0.1 m/sec. 

As Weinreb [1963] has shown, there is very little 
difference in digital processing time between the 
direct and autocorrelation routes to a power spec- 



trum when clipping is not permitted. What differ- 
ence there is, however, favors the autocorrelation 
function as a useful intermediary around which to 
reduce the fluctuation statistics of the signal by 
integration. 

One of the first to apply coherent digital processing 
to planetary echoes was Price [1959], In this first 
realization, 6-bit samples were taken of a filtered 
IF signal at 15 kc/s and stored on magnetic tape. 
Although permitting great flexibility in searching 
out a signal of unknown bandwidth and location in 
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Figure 4. Venus echo spectrum obtained by Millstone Radar 
of MIT at time of inferior conjunction, 10 April 1961. 
Analyzing slit width used was 0.2 c/s. 

The smooth curve represents an exponential approximation to the observa- 
tions. Transmitter frequency was 440 Mc/s. 

range and doppler frequency, this particular system 
required large quantities of computer time and would 
not be well matched to a continuing program. 

We return finally to an analog method which has 
proven useful in planetary radar. Brockman [1961] 
describes an application of a phase-locked receiver 
to the problem of measuring precisely the doppler 
shift of Venus echoes. The use of receivers which 
phase lock to an incoming signal and maintain long 
integration times in a feedback loop lias been com- 
monplace in telemetry applications for many years. 
As applied to the detection of weak signals from 
Venus (and Mercury) this method permits accurate 
determination of Doppler shifts of continuous wave 
echoes with a minimum of processing complexity. 
In practice it is limited to those cases where the 
returned Doppler spectrum is quite narrow, say a 
few tens of cycle per second at the most. Because 
it tends to follow the center of intensity of the spec- 
trum rather than the midpoint of the extremes, this 
method is probably less accurate than measurements 
made from a full spectrum. Nevertheless, in many 
applications its simplicity will make it attractive. 



3. Conclusion 

A number of spectral processing techniques have 
been described which range from fully analog to 
largely digital. Two avenues of approach to ob- 
taining a power spectrum of a weak or fluctuating 
signal are available, one direct and one via the auto- 
correlation function. The proper mix of these tech- 
niques in any given application will depend on the 
signal characteristics and on the objectives of the 
experiment. It seems clear, however, that as faster 
digital computers become available at reasonable 
prices, and as greater stability and more sophisticated 
objectives are sought, the digital techniques will be- 
come more commonly employed. 
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